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Abstract 

We employ the chiral nucleon-nucleon potential derived in ref.|l|| to study bound and scattering 
^ ' states in the two-nucleon system. At next-to-leading order, this potential is the sum of renor- 

P5 . malized one-pion and two-pion exchange and contact interactions. At next-to-next~to-leading 

order, we have additional chiral two-pion exchange with low-energy constants determined from 
pion-nucleon scattering. Alternatively, we consider the A(1232) as an explicit degree of freedom 
^\ ' in the effective field theory. The nine parameters related to the contact interactions can be de- 

termined by a fit to the np S- and P-waves and the mixing parameter ei for laboratory energies 
below 100 MeV. The predicted phase shifts and mixing parameters for higher energies and higher 
angular momenta are mostly well described for energies below 300 MeV. The S-waves are described 
as precisely as in modern phenomenological potentials. We find a good description of the deuteron 
properties. 
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1 Introduction 

Over the last years, effective field theory methods have been used to gain a better understanding of 
the two-nucleon interaction at low and intermediate energies. While at present these studies do not 
aim at substituting the highly successful "realistic" potentials build from meson exchanges (like the 
e.g. Bonn-Jiilich, Nijmegen, Argonne or RuhrPot potential), effective field theory (EFT) allows for 
a systematic and controlled expansion of observables in systems with two or more nucleons. Apart 
from dealing with the various scales appearing in nuclear systems, it is straightforward to implement 
the spontaneously and explicitely broken chiral symmetry of QCD as well as external probes in the 
EFT. The appearance of shallow bound states (or, equivalently, large scattering lengths) requires 
some method of resummation. There are essentially two ways of tackling this problem. One approach 
is to build the potential from EFT and employ it in a properly regularized Lippmann-Schwinger 
(or Schrodinger) equation. Alternatively, one can also do the expansion directly on the level of the 
scattering amplitude and resum the leading order, momentum-independent four-nucleon interaction. 
While at extremely low energies, much below the scale set by the pion mass, it is sufficient to consider 
four-nucleon interactions only (otherwise the rather successfull effective range expansion would not 
work), for typical nuclear momenta of the size of the pion mass, pions have to be included explicitely. 
While there has been much debate about the way how to treat the pions (perturbative versus non- 
perturbative pions), for the range of momenta to be considered here we believe that it is mandatory 
to include them at leading order and treat them nonperturbatively. The results discussed below give 
an a posteriori justification of this conjecture. In fact, our chiral potential at next-to-next -to-leading 
order in the power counting leads to results which are only slightly worse than the ones based on the 
so-called high accuracy modern potentials. 

In ref.||l| (referred to as I from here on) we constructed the two- and three-nucleon potential based 
on the most general chiral effective pion-nucleon Lagrangian using the method of unitary transforma- 
tions. For that, we developed a power counting scheme consistent with this projection formalism. In 
contrast to previous results obtained in old-fashioned time-ordered perturbation theory, the method 
employed leads to energy-independent potentials. This extends the power counting scheme origi- 
nally proposed by Weinberg [^ in a natural fashion. To leading order (LO), the potential consists 
of the undisputed one-pion exchange and two short-range four-fermion interactions. Corrections 
at next-to-leading order (NLO) stem from chiral two-pion exchange and additional contact terms 
with two derivatives. We extend the potential to next-to-next-to-leading order (NNLO) in two 
different ways. Apart from coupling constant and mass renormalization, one has only to consider 
the chiral two-pion-exchange potential (TPEP) with dimension two insertions from the pion-nucleon 
interaction. In contrast to ref.||3[, we take the novel low-energy constants (LECs) from systematic 
studies of pion-nucleon scattering in CHPT ^ ^, ^, |^.^ Alternatively, since most of these LECs 



are saturated by the A-resonance, one can also include the A explicitely.^ That only introduces 
the vrA^A coupling as a new parameter. It can be determined e.g. from an (spin-isospin) SU(4) 
relation to the pion-nucleon coupling. In our approach the exchanges of heavy mesons (like in the 
CD-Bonn or Nijm93 potential) or parametrizations thereof (like in the AV18 or NijmI,II potential) 
show up as contact terms, but we refrain from a comparison of similarities and differences to the var- 
ious potential model approaches at this point. A more detailed introduction of our method is given 
in I. Here, we will be concerned with numerical results obtained on the basis of that two-nucleon 



*^We already remark here that the determination of these parameters from fitting the invariant amplitudes inside the 
Mandelstam triangle is favored by our fits. 

*^ Strictly speaking, one should keep also the dimension two -kN operators and subtract the A contribution. Since an 
explicit A is, however, dynamically different from integrating it out, we refrain from doing this. The uncertainty due to 
this procedure is small in most partial waves. 



potential at NLO and NNLO. For that, we first have to renormahze the potential and discuss the 
appropriate cut-off regularization. Having done that, it is straightforward to solve the corresponding 
Lippmann-Schwinger equation for bound and scattering states. In total, we have nine parameters 
related to the four-fermion contact terms. These parameters can be uniquely fixed from fitting np 
partial waves at low energies, more precisely the two S-, four P-waves and the mixing parameter ei. 
Predictions for higher energies and the higher partial waves and deuteron properties arise (we refrain 
from fine tuning some of the parameters in the mixed ^Si — ^ Di waves to get the exact deuteron 
binding energy as it is done in conventional boson-exchange models). 

To put our calculations in a better perspective, we briefly review the status of previous works (directly 
related to the results presented) without going into details. The idea of constructing the potential 
from effective field theory was put forward by Weinberg |Q] in the early nineties. This was taken up by 
van Kolck and collaborators and culminated in ref.jy]. In that paper, the two-nucleon potential was 
constructed at next-next -to-leading order ^ and analyzed. It contained one- and two-pion exchange 



graphs (based on time-ordered perturbation theory diagrams) and a host of contact interactions. An 
exponential cut-ofF at each vertex was introduced to tame the high energy behaviour. Fierz reordering 
was not used so that global fits with 26 parameters had to be performed and first numerical results 
for phase shifts and deuteron properties were obtained. While promising, these results were not as 
precise as the ones usually obtained in boson-exchange models. Park et al. M considered a series of 
interesting applications mostly related to the deuteron properties. They restricted themselves to one- 
pion-exchange and contact interactions for the relevant phases ^Sq,'^ S*!,^ Di and the mixing parameter 
ei. Peripheral nucleon-nucleon scattering was considered by the Munich group Q, including two- 
pion exchange graphs with insertions from the pion-nucleon Lagrangian of dimension two based on 
dimensionally regularized Feynman graphs. This lead to an improved description of some higher 
partial waves (for similar results, see the work of |lO|). Since the potential was only considered 
perturbatively, the bound-state problem could not be addressed. A different counting scheme was 



proposed by Kaplan, Savage and Wise |11|. In that framework, the low partial waves and mixing 
parameters were considered to NLO and NNLO as well as many deuteron properties. A NNLO 
calculation of the ^Si —^ Di transition potential matrix element was recently presented |12]. Further 



NNLO investigations for other partial waves in the framework of the KSW approach are under way 
by the Seattle group, see ref.in 



Our manuscript is organized as follows. In sec.|2|, we explicitely give the renormalized potential at NLO 
and NNLO. This follows directly from the potential derived in I.^ We also discuss the regularization 



procedure necessecary to render the (iterated) potential finite. The Lippmann-Schwinger equation 
underlying the calculation of the scattering and bound states is briefly discussed in sec.|3|. The fitting 
procedure to determine the LECs and the accuracy of the fits are detailed in secQ. This involves the 
low partial waves at (kinetic) energies (in the lab frame) below 100 MeV. Results for the low partial 
waves at higher energies, for the higher partial waves and the deuteron (bound state) properties are 
displayed and discussed in sec.||. Our findings are summarized in sec]6|. The appendices contain 
details on the renormalization procedure, the inclusion of the A-resonance and a collection of other 
useful formulas. 



*'^The leading effects of the A (1232) resonance were also incorporated in that work. 
**Note that the NNLO potential was not given in I. 



2 The renormalized potential 

In I, we derived the two-nucleon potential at next-to-leading order. To leading order, it consists 
of the one-pion-exchange potential (which is called the OPEP) and two S-wave (non-derivative) 
contact interactions. The latter are parametrized by the coupling (low-energy) constants Cs and Ct- 
The OPEP constitutes the longest range part of the nucleon-nucleon interaction and it completely 
dominates the high partial waves. However, the tensor potential related to one-pion-exchange be- 
comes unrealistic at short distances. Chiral symmetry enforces a (pseudovector) derivative coupling 
in harmony with the Goldstone theorem. Note also that while the four-fermion interactions are 
pointlike on the level of the Lagrangian, they get smeared out by the regularization procedure to be 
discussed below. At next-to-leading order, we have three distinct contributions. The first ones are 
the one-loop self-energy and vertex corrections to the OPEP, cf. figs. 3, 4 in I. Similarly, it is manda- 
tory to also work out the one-loop corrections to the contact terms ~ Cs, Ct, cf. fig. 6 in I. Thirdly, 
there are the genuine two-pion exchange diagrams, leading to the so-called TPEP, cf. fig. 2 in I. This 
TPEP is the second longest range component of the A^A'^ interaction. Here, chiral symmetry leads 
to the NNtttt vertex which in turn leads to the so-called triangle and football diagrams. These are 
usually not accounted for in boson-exchange models. One notices immediately that some of the loop 
corrections contain UV divergences. Since we are dealing with an effective field theory, we can use the 
standard order-by-order renormalization machinery for the Lagrangian, which was first developed in 
the context of chiral perturbation theory. All these divergent contributions can be renormalized in 
terms of three divergent (momentum space) loop functions, 

Jo = / — , J2=/ Idl, Ja= l^dl. (2.1) 

Je I Jo JO 

Clearly J2 ( J4) is quadratically (quartically) divergent while Jq diverges logarithmically for large 
momenta and e is an IR regulator. Note that we could introduce other forms of the divergent integrals, 
which depend on one dimensionful scale (denoted here by e). Another choice of defining this divergent 
loop integrals would give different finite subtractions but leads to the same non-polynomial terms 
in the potential. The precise renormalization procedure is discussed in app.^. Notice also that 
different to standard use, we perform the renormalization of the potential before regularizing the LS- 
equation as discussed below. Such an additional regularization is needed since loop as well as contact 
term contributions grow quadratically in momenta. More precisely, at NNLO some of the TPEP 
contributions even grow with the third power of momentum. Such an UV behaviour of the potential 
is, of course, unacceptable and requires additional regularization. The solutions of the LS-equation 
are adjusted to observables. Therefore, the renormalization constants are fixed and related to the 
regularization procedure in the LS-equation. This procedure is justified since the potential by itself 
is not an observable but rather the bound and scattering states obtained from the LS-equation. This 
should always be kept in mind. 

We consider first the contact terms of the two-nucleon potential. To the accuracy we are working, 
the matrix-element of the potential in the center-of-mass system (cms) for initial and final nucleon 
momenta p and p', respectively, takes the form (note that due to the choice of the cms, a reduction 
in the number of terms follows as has been shown in I) 

y(contact) _ y(0) _|_ y(2) 

1/(0) = Cs + CTai-a2, 

F(2) = Ci f + C2P + (C3 f + C4 P)(ai • <T2) + iCs ^ (<Ti + a2) ■ {g x k) 

+ Cq {q- ai){q- 02) + C7 (fc • ai)(fc • 02) , (2.2) 



with q = p ' — p and k = {p + p ')/2- These terms feed mto the matrix-elements of the two S-waves 
(^S'o,'^ Si), the four P-waves i^Pi,^ Pi,^ P2,^ Po) and the lowest transition potential i^Di —^ Si) in the 
following way: 



V{^So) = 47r {Cs - 3C7t) + ^ (4C7i + C2 - I2C3 - 3C4 - 4^ - C^) (/ + p>^ 



C'iSo + C^iSo(p'+p"), (2-3) 



y(35i) = 47r (Cs + Ct) + I (12C7i + ZC2 + 12C^ + 3C4 + ACq + C7){p^ + p'^) , 

= C^s,+C^sAp'+p")^ (2-4) 

2tt 
Vi^Pi) = — i-iCi + C2 + I2C3 - 3C4 + 4C76 - C7) ipp') = Cip, (ppO , (2.5) 

O-jr 

Vi^Pi) = — (-4Ci + C2- 4C3 + C4 + 2C5 + 4C76 + C7) (ppO = Csp^ ipp') , 

(2.6) 

v^('P2) = Y (-4Ci + C2- 4C3 + C4 + 2C5) (ppO = ap, (ppO , (2.7) 

2tt 
y(3po) = Y (-4Ci + C2- 4C3 + C74 + 4C5 + 12C6 - 3C77) (pp') = Cip, {pp') , 

(2.8) 
y(3Di-3 5i) = '^{iC, + Cr)p'' = an,^ss,p'\ (2-9) 

n'^i-^Di) = ^(4C6 + C7)p2 = C3B,_3s,P^ (2.10) 

with p = \p\ and p' = |p'|. These nine constants are not fixed by chiral symmetry but can be 
determined by a fit to these lowest partial waves and one mixing parameter as detailed below. We 
have already given the appropriate partial wave decomposition for the low-energy constants (LECs) 
Cs, Ct, Ci, . . . , C7 here. From each of the two S-waves, we can determine two parameters, whereas 
the four P-waves and the mixing parameter ei contain one free parameter each. Of course, we have 
to account for the channel coupling in the mixed triplet partial waves. It is also important to note 
that once the C2S+1;. have been determined, the original Cs, Ct, Ci, . . . , C7 are fixed uniquely. We 
remark that the values for the C23+11. are renormalized quantities, see app.^. 

Consider now OPEP and TPEP. After vertex and coupling constant renormalization, as detailed in 
app.^, we find the following expressions in terms of the renormalized quantities at NLO: 

T.OPEP f9AY^ ai-qa2-q . 



^NL^o" = -^^.L{q){^Ml{f>g\-Ag\-l)+q\2^g\-mg\-l) + ^^_ 



4 A/f4 
2 



^9\ 



647r2/^ 



with 



L{q) lai ■ qd2 ■ q- q^ Bi ■ 5^ + P(fe, q) , (2.12) 

L(,) = ivWT7ln^^?^|^^, (2.13) 

q * 2iw^ 

and we have set q = \q\. Furthermore, /vr = 93MeV is the pion decay constant, M^^ = 138.03 MeV 
the pion mass and gA = 1-26 for the axial-vector coupling.^ Furthermore, P{k,q) is a polynom 

* Note that we have changed our conventions as compared to I for the pion decay constant, the isospin generators 
and the relative momentum q. 



in momenta of at most second degree and has the same structure as the expressions in eq.(p.2|). 
More precisely, after performing the partial decomposition, P{k, q) leads in each partial wave to 
polynoms in p and p' of at most second degree. Thus, its explicit form is of no relevance here, since 
it only contributes to the renormalization of the couplings Cs,CT,Ci, as detailed in app.^. The 
TPEP agrees with the one given by the Munich group. It is important to stress the differences to 
the calculation of ref.Q. While there the potential was treated perturbatively, we iterate it to all 
orders in a Lippmann-Schwinger (LS) equation. Second, we use a cut-off regularization within the 
LS-equation and not dimensional regularization on the level of the diagrams. Of course, for the 
peripheral partial waves, the iteration is not of importance. We are, however, more ambitious in that 
we want to get a description of all partial waves as well as of the bound state (deuteron) properties. 
The TPEP at NNLO has also been given in ref.|p using dimensional regularization. Within our 
renormalization scheme, it reads j^n 

9 a ,^ ^^r 35iM5 



1287rm/4 



(n • r,) {-^£^, + ml + 2,2 - g\ml + 3q'))i2Ml + q')Aig)^ 



+ ^^l^f, ((^1 • ^K'?^ • 9 ) - q\^i ■ ^2)) {2Ml + q^)A{q) 
^ :(n • T2) (<?i • g)(<T2 • q) - q^fyffi ■ 0^2) 



327r/4 



: |(C4 + ^)ml + q^) - ^{WMl + Sg^)} A{q) 



^9 A i^ff^+ff^).^p',,f)^2Ml + q')A{q) 



647rm/4 



647rm/^ 



with 



(n • fa) i (<Ti +a2)-{p'xp) {AMI + Q)Mq) + P{K q) , (2-14) 



^('^) = ^^^^^^^2l^- ^'-''^ 



The polynom P{k, q) plays a similar role as P{k, q) in eq. (|2.12| ). Note that we have used the standard 
dimension two vriV Lagrangian as in ref.||9[. For the LECs 01^3^4 we should take the values obtained from 
fitting vrA^ phases in the threshold region, see e.g. ref.Q or, alternatively, from fitting the invariant 
amplitudes inside the Mandelstam triangle, i.e. in the unphysical region Ij. The so determined 
parameters are only slightly different, but these small differences will play a role later on. For 
example, the LECs 01^3^4 from fit 1 of ref.|^] are ci = — 1.23GeV~^, C3 = — 5.94GeV~^ and C4 = 
3.47 GeV~^. A recent investigation of the subthreshold amplitudes |Q] leads to slightly different 
values, ci = —0.81 GeV~^, C3 = — 4.70GeV~^ and C4 = 3.40 GeV~^. It is this latter set we will use in 
the following. These values are also consistent with the recent determination from the proton-proton 
interaction based on the chiral two-pion exchange potential [^. 

The total renormalized potential at NLO is now given as the sum of the OPEP, TPEP and contact 
potentials as given in eqs.( 2.ll]j2.12 ,2^). At NNLO, we have to add the additional TPEP from 



Note that in I we have used the power cou nting such that Q/m ~ Q /A^, where Q corresponds to the low momentum 



scale. Accordingly, the 1/m corrections in eq.(2.14) are smaller than those given by the ci, C3, and C4 terms and therefore 



contribute to the N'^LO potential. In fact, this follows also if one compares the numerical vaues of these two types of 
corrections (1/m versus a). We have decided here to keep these 1/m corrections, since otherwise one cannot directly 
use the values of the LECs ci_3,4 as determined from the ttN sector in the presence of the 1/m terms. 



eg. ( ^.14 ). Alternatively, we have also included the leading contribution from the A(1232) resonance. 



This is formally a NNLO contribution as demanded by the decoupling theorem. The pertinent 
equations are collected in app. ^. We only remark that there is one new parameter, namely the AA^vr 
coupling constant. It can be determined from the SU(4) relation 

3 gAm 
9AN. = ^^-^ (2-16) 

making use of the Goldberger-Treiman relation (GTR) (higher order corrections to the GTR are not 
considered here). This leads to g/^N-n = 27.2, as detailed in app. p. The corresponding theory is 
called NNLO-A. The chiral TPEP has a different high momentum behaviour as in the strict NNLO 
approach because of the explicit A-propagator (assuming that we only include the A to generate 
the NNLO TPEP and no additional ttN contact interactions). Only for small momenta NNLO and 
NNLO~A are essentially the same (resonance saturation). At this point, we already mention that 
the NNLO TPEP is too strong at short distances (large momenta). To the order we are working, we 
have counterterms in the S- and P-waves to balance this, but e.g. not in the D-waves. This will be 
discussed in more details below. 

Since the potential is only meaningful for momenta below a certain scale, it needs regularization. 
Stated differently, the large momentum behaviour of the potential is not correct. To the order we 
are working, the contact interactions diverge quadratically for large momenta and some components 
of the TPEP grow even with the third power of momentum. As it is appropriate in effective field 
theory, we regularize the potential. That is done in the following way: 

V{p, p')^ fnip) Vip, p')fR{p'), (2.17) 

where fR{p) is a regulator function chosen in harmony with the underlying symmetries. In what 
follows, we work with two different regulator functions, 

fp^^ip) = eiA^-p'), (2.18) 

/r°°(p) = exp(V7A'"), (2.19) 



with n = 2,3,.... The sharp cut-off is most appropriate to the projection formalism. For the 
calculation of some observables, however, it cannot be used since at p,p' = A it leads to discontinous 
derivatives. For very large integers n the exponential cut-off approximates the sharp one. Throughout, 
we work with n = 2. To the order we are working, the choice n = 1 has to be excluded since the 
contact terms of order p'^,p' would be modified. For n = 2, the error we make is beyond the accuracy 
of the order we are calculating. We are now in the position to calculate observables with this potential. 

3 Bound and scattering state equation 

In ref.|^], the Schrodinger equation was solved after the chiral NN potential had been transformed from 
momentum into co-ordinate space. We consider it more natural to work directly in momentum space. 
The corresponding equation describing the bound and scattering states is the Lippmann-Schwinger 
equation. Here, we briefly discuss it to keep the manuscript self-contained. For a more detailed 
exposition concerning also methods of solving the LS-equation, we refer to the monograph |15|| . 

The LS-equation (for the T-matrix) projected into states with orbital angular momentum /, total 
spin s and total angular momentum j is 



oo 



dp P .rsj r I ii\ 1 



mp\p) = ^iW^p) + E X w ^"^"(^''^")^^37^3-^"/^'» ' (3.1) 



with T] -^ 0+. In the uncoupled case, I is conserved. The partial wave projected potential Vif\{p',p) 
is obtained as follows. We first rewrite the potential V in the form ||l^ 



V = Vc + Va ai-a2 + VsL i i^{^i + ^2) ■ {k x q) + V^l <?! • (g x k) a2-{qx k) 
+ Vaq (cti -q) {^2-q) + Vak (<?i • k) (ct2 • k) 



(3.2) 



with six functions Vc{p,p' , z), . . . , Va-kip,p' , z) depending on p = \p\, p' = \p'\ and the cosine of the 
angle between the two momenta z. These functions may depend on the isospin matrices r as well. 
For j > 0, one obtains in the usual Isj representation 

(jOJlFlJOi) = 2t, j dz [Vc - W^+p'^p'iz' - l)V^L - q^Va<t - k'V^k} P,{z) , 



{nj\V\jlj) = 27tJ dz 



.'2^2 



Vc + V^ + 2p'pzVsL - p' V(l + 3z')Kl + ^k'V^, + -rq'V^k 



xP,{z) + 



.'2^2 



p'p VsL + 2p''p'zV^L - 2p'p (y,, - -Kfc 



x(P,_i(z) + P,+i(z)) , 



(j ± 1, lj\V\j ± 1, Ij) = 2-K |_^ dz |pV 

X Pj[z) + 



-V5L± 



VV^VaL + Vuq - -^ak 



2j + l 

Vc + Va+ V'PZVSL + P'V(1 - Z^)V^L 



1 



"2j + l 



.'2 2 



1, 



2p'VKl - {p' +p')iy„q + T^afc 



Pmi^) 



{j ± 1, lj|y I j T 1, Ij) = ^^|I^ll27r 1'^ dz I - p'p (4K, - Kfc)P,(^) 



+ 



+ 



2p p 



2i + 
2p'^p2 



2j + 



^VaL+p'\2V„q + ]^V„t 
' 2/ 1 



PlTliz) 
P3±li.A 



(3.3) 



Here, -Pj(z) are the conventional Legendre polynomials. For j = the two nonvanishing matrix 
elements are 

(000|y|000) = 2tt j dz[Vc-W„+p'^p'{z^-l)V„L-q^V„q-k'V„k} , 

(110|y|110) = 27ry dz{zVc + zV,+p'p{z^-l)VsL+p'^p'z{l-z^)V„L (3.4) 

- ((P'' + P')z - 2p'p) V^g - ^ [{p'^ + p^)z + 2p'p) V„ 



Note that sometimes another notation is used in which an additional overall "-" sign enters the 
expressions for the off-diagonal matrix elements with / = j + 1, I' = j — 1 and I = j — 1, I' = j + 1. 

The relation between the on-shell S- and T-matrices is given by 



Snip) = Si'i - -^P^T^'lip) 



(3.5) 



where p denotes the two-nucleon center-of-mass three-momentum. The phase shifts in the uncoupled 
cases can be obtained from the 5-matrix via 

5°^ = exp(2.^^), 5g = exp (2.^^) , (3.6) 

where we have used the notation 5p. Throughout, we use the so-called Stapp parametrization |]l 
of the S'-matrix in the coupled channels {j > 0) : 

g ^ ( ^l-ij-i ^i-^J+^\ = ( cos(2e)exp(2i5]li)^ i sin {2e) e^p {i5]l^ + iS^/_^^) 



S^+ij-i 'S'jii j+i / V ^ sin (2e) exp (idji^ + i6]{^) cos (2e) exp (2i(5]:^^) 

The bound state is obtained from the homogeneous part of eq.(0) and obeys 



(3.7) 



Ed - pym y" Jo {2TTy '•' 

with s = j = 1 and / = /' = 0, 2 and E^ denotes the bound-state energy. The LS-equations for 
the scattering and the bound state(s) are solved by standard Gauss-Legendre quadrature. It goes 
without saying that in this discussion, the potential is to be understood in its regularized form as 
detailed in eqs.^T^, |2T|, |1|). 



4 The fits 

In this section we discuss the determination of the various coupling constants. The leading OPEP is 
constrained by chiral symmetry, its strength is given by the pion-nucleon coupling constant, which 
due to the Goldberger-Treiman relation is proportional to qa/ f-rr- For the NNLO-A approach, we in 
addition have the strong pion-nucleon- A coupling constant. 

To pin down the nine parameters C5, Ct, Ci, . . . , C7 we do not perform global fits as done in ref.Q. 
Rather we introduce the independent new parameters as given in eqs.(2.3- 2.10| ). So to leading order. 



the two S- waves are depending on one parameter each. At NLO, we have one additional parameter 
for ^Sq and ^Si as well as one parameters in each of the four P-waves and in ei. We thus can fit each 
partial wave separately, which makes the fitting procedure not only extremely simple but also unique. 
Of course, in case of the triplet coupled waves i^Si —^ Di,^ P2 —^ -F2) the fitting is performed for the 
corresponding 2x2 S-matrix parametrizations (parametrized by two partial waves and one mixing 
parameter). At NNLO (NNLO-A), we have no new parameters, but must refit the various contact 
interactions due to the TPEP contribution in all partial waves. We have used two different methods to 
fix the LECs of the contact interactions. First, we fit to the phase shifts of the Nijmegen partial wave 
analysis |2^ for laboratory energies smaller than (50) 100 MeV at (NLO) NNLO. P^ Alternatively, 
we used the effective range parameters for the phases ^5*0, ^Si and ei to fix Ci^^, (715^,(73 s^, C3 5^, 
and Cs£)^_ssi- I^^ what follows we will mark the corresponding potentials by a "•" if the LECs were 
fixed from the effective range parameters. We have found that the observables (phase shifts) depend 
rather weekly on what fitting procedure is used. The values of the LECs do not change much except 
for some C"s at NLO, where significant variations were observed. Note further that in the S-waves, 
especially in ^Sq, isospin breaking effects like e.g. the charged to neutral pion mass difference, are 
known to be important. We do not consider such effects in this work and thus take an average value 

Note that equally well we could fit directly to the data. However, for easier comparison with other EFT calculations, 
wc use the Nijmegen PSA to simulate the data. 



for the pion mass. The actual values of the S-wave LECs seem to be rather sensitive to the choice 
of the pion mass. For the phases ^5*0, ^5*1 and ei we use the errors as given in ref.[^, for all other 
partial waves we assign an absolute error of 3%. This number is arbitrary, but taking any other value 
would not change the fit results, only the total x^. To perform the fits, we have to specify the value of 
the cut-off A in the regulator functions as defined in eqs.( 2.18 . 2.1S| ). At NLO, for any choice between 
380 MeV and 600 MeV, we get very similar fits (a very shallow x^ distribution in each partial wave). 
At NNLO, this shallow distribution turns into a plateau, which shifts to higher values of the cut-off. 
We can now use values between 600 MeV and 1 GeV. The results using the sharp or the exponential 
cut-off are very similar. For illustration, we consider the sharp cut-off with A = 500 MeV at NLO 
and A = 875 MeV at NNLO. We show in fig. || the absolute quadratic deviations of the fit to the 
Nijmegen phase shift analysis (PSA), defined by ((5^* — 6^^ )'^. Consider first the S-waves. Both in 
^Sq and ^Si, we observe a clear improvement when going from LO to NLO to NNLO (NNLO-A). 
A similar pattern holds for the P-waves and ei, although the differences between NLO and NNLO 
are somewhat less pronounced. Note that the ^Pq wave is very senstive to the value of the pion 
mass, therefore the slightly better NLO fit should not be considered problematic. The corresponding 
parameters of the coupling constants are collected in table |i|. 





C'l^o 


CiSo 


Cisi 


Cssi 


CiPi 


Cip^ 


Cap, 


Csp^ 


CsDi-ss-^ 


NLO 


-0.134 


1.822 


-0.130 


-0.393 


0.344 


-0.394 


1.335 


-0.1907 


-0.0317 


NLO* 


-0.0928 


2.125 


-0.102 


0.0243 


0.344 


-0.394 


1.335 


-0.191 


-0.0357 


NNLO 


-4.249 


11.945 


-6.508 


11.293 


-2.045 


-7.061 


-2.832 


-8.056 


-3.424 


NNLO* 


-4.246 


11.943 


-6.655 


11.279 


-2.045 


-7.061 


-2.832 


-8.056 


-3.516 


NNLO-A 


-5.731 


13.823 


-0.977 


3.446 


-2.494 


-8.188 


-2.993 


-8.431 


-0.562 



Table 1: The values of the LECs as determined from the low partial waves. We use a sharp cut-off 
with A = 500 MeV and 875 MeV at NLO and NNLO, respectively. The Ci are in lO'^ GeV'^ while the 
others are in 10^ GeV~^. The parameters of the NLO, NNLO and NNLO-A potentials are obtained 
from fitting to the Nijmegen PSA. The LECs of the NLO*, NNLO* potential in the ^^o and ^Si -^ Di 
channels are fixed to reproduce exactly the effective range parameters. 

To illustrate the dependence on the cut-off, we show in fig.0 the running of the two (three) couplings 
in the ^Sq {"^Si) channels at NNLO (note that the third parameter in ^Si comes in via the mixing 
with the ^Z?! wave). We notice that the variation of these LECs over a wide range of cut-offs is rather 
modest. We also mention that using the vrA^ parameters from refs.Q, [s], g leads to a considerably 
worse x^ iii the fits. We take that as an indication that the determination of the q based on the 
method employed in ref.[0] is more reliable than fitting to vrA^ phase shifts (as long as one works to 
third order in the chiral expansion). We remark that using the parameters of ref.[0], the deuteron 
binding energy Eij_ comes out as 



NLO : ^rf = -2.175MeV 
NNLO : ^rf = -2.208MeV 



(4.1) 



i.e. the NNLO result is already within 7.5 permille of the experimental number. Fine tuning in 
the parameters in the deuteron channel would allow to get the binding energy at the exact value of 
—2.224575(9) MeV without leading to any noticeable change in the phase shifts. We later consider 
the deuteron channel separately with an exponential regulator. This will lead to an improved binding 
energy but no attempt is made to match the exact value in all digits. 
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At this point we would like to comment on the increase in the cut-off values when going from NLO to 
NNLO (NNLO-A). Consider first the leading order result. Lepage [23| has pointed out that inclusion 



of the one-pion exchange does not lead to a remarkable increase in the cut-off values compared to a 
pionless theory. In order to fit the phase shifts in the S-waves one should choose the cut-off below 
500-600 MeV even if the contact interactions with two derivatives are taken into account (within 
our power counting scheme such contact interactions contribute first at next-to-leading order and 
are of the same size as the lading two-pion exchange terms). Lepage assumed that such a low value 
of the cut-off is due to the missed physics assosiated with the two-pion exchange. So, naively, one 
would expect that the inclusion of the leading two-pion exchange contributions at NLO would allow 



to take larger cut-off values. However, that does not happen. This was also pointed out in ref.|24|. 
According to our analysis only at NNLO, after the subleading two-pion exchange contributions are 
taken into account, one can increase the cut-off up to 800 to 1000 MeV. The inclusion of the dimension 
two operators of the pion-nucleon interaction at NNLO encodes some information about heavy meson 
exchange as well as virtual isobar excitations, as discussed in detail in ref.Q. In this work we were able 
to separate the leading effects of A-resonance (NNLO-A). The clear increase in the cut-off values 
when going from NLO to NNLO-A indicates the importance of physics assosiated with heavier mass 
states like e.g. the A-resonance. Our NNLO (NNLO-A) TPEP is sensitive to momentum scales 
sizeably larger than twice the pion mass (as it would be the case for uncorrelated TPE) and delta- 
nucleon mass splitting. Consequently, the cut-off has to be chosen safely above these scales, say above 
500 MeV (with the sharp regulator). The upper limit of about 1 GeV is related to the cancellations 
in the S-waves (fine-tuning), for too large values of A it is not longer possible to keep this intricate 
balance. It is, however, comforting to see that including more physics in the potential leads indeed 
to a wider range of applicability of the EFT. 

Finally, we need to discuss one further topic. Performing the fits, we have found two minima in both 
the ^5*0 and the "^Si channel. This is not unexpected and can easily be understood in the case of a 
pionless theory. In that case, exact analytical calculations are possible. As was shown in ref.|25], the 



requirement of reproducing exactly the scattering length and the effective range leads to the following 
conditions for the LECs in the ^5*0 partial wave. 



mr_ ^ f m \^ C.s,{2 - C.s.h) 

327r2 \lQTT^a ^) (Cisjs-l)^ ' (2^)3A ' 



+ T^TZvTT ' (4-2) 



,^dq„^ 771 A" , , 

In = -m TTT^g"-' = -t;^^— • (4.3) 



where 

~~^ lo (2^'' ~ " (2^ V ■ 

The second equation in ( |4.2| ) can be further simplified if a ^ 1/A. Then, one obtains the following 
quadratic equation for Cig^, 

m { 4\ Cig.,(2-Cig,/3) 



which leads to 

Note that the existence of real solutions for Cig^ requires that A < 16/(7rr) ~ 376 MeV. Such a 
situation with two solutions appears also in the NLO and NNLO theory with pions. At NLO, we find 
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very similar predictions for the phase shifts and observables as well as a very similar quality of the 
fits in the ^Sq and '^Si —^ Di channels for both solutions, see the upper panel in fig.|3|. So there is no 
real criterion to prefer one of these solutions. However, at NNLO, the behaviour of the phase shifts 
at higher energies differs quite remarkbaly as it is illustrated in the lower panel of fig.^. Also, the x^ 
for these two solutions differs typically by factors of 2. ..10. In what follows, we will only discuss the 
best solution at NNLO. 

5 Results and discussion 

Having determined the parameters, we can now predict the S- and P-waves for energies larger than 
100 MeV and all other partial waves are parameter free predictions for all energies considered. Since 
at the laboratory energy of 280 MeV the first inelastic channel opens, we only plot the phase shifts 
up to E'lab = 300 MeV. We show the results using the sharp regulator function with A = 500 and 
875 MeV for NLO and NNLO, as discussed in sec.^. This is a best global fit for all partial waves. 
Using the exponential cut-off, the curves come out very similar. In fact, for the calculation of some 
deuteron observables, we have to use the exponential regulator. Obviously, the corresponding cut-offs 
are somewhat larger than in case of the sharp regulator (to preserve the area to be integrated). For 
the NNLO-A approach, we only show some selected partial waves in a separate paragraph. 

5.1 Phase shifts 

5.1.1 S waves 

In fig.H we show the two S-waves at LO, NLO and NNLO for the cut-off values given above. Clearly, 
the lowest order OPEP plus non-derivative contact terms is insufficient to describe the ^Sq phase (as 
it is well-known from effective range theory and previous studies in EFT approaches). The much 
more smooth ^Si phase is already fairly well described at leading order. For energies above 100 MeV, 
the improvement by going from NLO to NNLO is clearly visible. The corresponding values of the 
S-wave phase shifts at certain energies are given in tables |2|,y. For comparison, we also give the 
results of the Nijmegen and VPI PSA [21] and of three modern potentials (Nijmegen 93 |22], Argonne 



V18 \2t] and CD-Bonn [^]). Our NNLO result for ^Sq is visibly better than the one obtained in 



ref.|l|]. 



It is also of interest to consider the scattering lengths and effective range parameters. The effective 
range expansion takes the form (written here for a genuine partial wave) 

p cot((5) = — + - r/ + V2P^ + V3p^ + Vip^ + 0{p^°) , (5.1) 

a 2 

with p the nucleon cms momentum, a the scattering length and r the effective range. It has been 
stressed in ref. p^] that the shape parameters are a good testing ground for the range of applicability 
of the underlying EFT since a fit to say the scattering length and the effective range at NLO leads to 
predictions for the Vi . ^^ In table ^, we present our results for the S-waves in comparison to the ones 



obtained from the Nijmegen PSA. Note that in the coupled channel we have used the so-called Blatt 
and Biedenharn parametrization of the S-matrix in order to be able to compare our findings with 



those of the ref. 1 19]. To show that our results are stable and do not depend on the fitting procedure 



we present the NNLO predictions obtained with two different methods of fixing the LECs as discussed 



*^ We remark in passing that the so-called "low-energy theorems" discussed in that paper do not qualify under this 
title. For a pedagogical discussion on this point, see e.g. ref. ]29(]. 
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above. Overall, the agreement is quite satisfactory. The same holds for the NLO results, but as stated 
before, these are more sensitive to the fitting procedure. 



i^iab [MeV] 


NNLO* 


NNLO 


Nijm PSA 


VPI PSA 


Nijm93 


AV18 


CD-Bonn 


1* 


147.735 


147.727 


147.747 


147.781 


147.768 


147.749 


147.748 


2 


136.447 


136.450 


136.463 


136.488 


136.495 


136.465 


136.463 


3 


128.763 


128.781 


128.784 


128.788 


128.826 


128.786 


128.783 


5* 


118.150 


118.196 


118.178 


118.129 


118.240 


118.182 


118.175 


10* 


102.56 


102.67 


102.61 


102.41 


102.72 


102.62 


102.60 


20 


85.99 


86.21 


86.12 


85.67 


86.35 


86.16 


86.09 


30 


75.84 


76.14 


76.06 


75.46 


76.40 


76.12 


75.99 


50* 


62.35 


62.79 


62.77 


62.12 


63.36 


62.89 


62.63 


100* 


42.33 


43.06 


43.23 


42.98 


44.33 


43.18 


42.93 


200 


19.54 


20.68 


21.22 


20.88 


22.82 


21.31 


20.88 


300 


4.15 


5.58 


6.60 


5.08 


8.44 


7.55 


6.70 



Table 2: ^5*1 np phase shift for the global fit at NNLO (sharp cut-off, A = 875 MeV) compared 
to phase shift analyses and modern potentials. The parameters of the NNLO potential are fixed 
by fitting the Nijmegen PSA at six energies (-Biab = 1, 5, 10, 25, 50(, 100) MeV). These energies are 
marked by the star. The parameters of the NNLO* potential are choosen to reproduce exactly the 
scattering length and the effective range as described in the text. 



i^lab [McV] 


NNLO* 


NNLO 


Nijm PSA 


VPI PSA 


Nijm93 


AV18 


CD-Bonn 


1* 


62.071 


62.063 


62.069 


62.156 


62.065 


62.015 


62.078 


2 


64.472 


64.469 


64.573 


64.573 


64.460 


64.388 


64.478 


3 


64.671 


64.671 


64.762 


64.762 


64.650 


64.560 


64.671 


5* 


63.659 


63.663 


63.708 


63.708 


63.619 


63.503 


63.645 


10* 


60.02 


60.03 


59.96 


60.00 


59.94 


59.78 


59.97 


20 


53.66 


53.68 


53.57 


53.77 


53.54 


53.31 


53.56 


30 


48.55 


48.58 


48.49 


49.00 


48.42 


48.16 


48.43 


50* 


40.49 


40.54 


40.54 


41.66 


40.38 


40.09 


40.37 


100* 


26.30 


26.38 


26.78 


27.86 


26.17 


26.02 


26.26 


200 


7.63 


7.76 


8.94 


7.86 


7.07 


8.00 


8.14 


300 


-6.41 


-6.24 


-4.46 


-5.55 


-7.18 


-4.54 


-4.45 



Table 3: ^5*0 np phase shift for the best fit at NNLO (sharp cut-off, A = 875 MeV) compared to phase 
shift analyses and modern potentials. The parameters of the NNLO potential are fixed by fitting the 
Nijmegen PSA at six energies (.Biab = 1, 5, 10, 25, 50, 100 MeV). These energies are marked by the 
star. The parameters of the NNLO* potential are choosen to reproduce exactly the scattering length 
and the effective range as described in the text. 
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We also remark that the momentum expansion for ei discussed in ref . |28f| depends sensitively on the 
parametrization one uses for the coupled triplet waves and we thus refrain from discussing this issue 
any further. ^^ 





a [fm] 


r[fm] 


V2 [fm3] 


^3 [fm^] 


Vi [fm^] 


^So NNLO* 


-23.739 


2.68 


-0.61 


5.1 


-29.7 


1^0 NNLO 


-23.722 


2.68 


-0.61 


5.1 


-29.8 


^So NPSA 


-23.739 


2.68 


-0.48 


4.0 


-20.0 


35i NNLO* 


5.420 


1.753 


0.06 


0.66 


-3.8 


35i NNLO 


5.424 


1.741 


0.05 


0.67 


-3.9 


35i NPSA 


5.420 


1.753 


0.04 


0.67 


-4.0 



Table 4: Scattering lengths and range parameters for the S~waves at NNLO (global fit) compared 
to the Nijmegen PSA (NPSA). The values for ^2,3,4 in the ^^o channel are based on the np Nijm II 



potential and the values of the scattering length and the effective range are taken from the ref. [18]. 
The effective range parameters for the "^Si —^ Di channel are discussed in [19|. 



For completeness, we collect the experimental values for the S-wave scattering lengths and efffective 
ranges: 



^^0 



at 



(-23.758 ±0.010) fm, 
(5.424 ± 0.004) fm, 



n 



= (2.75 ± 0.05) fm , 
(1.759 ± 0.005) fm 



(5.2) 
(5.3) 



5.1.2 P waves 

In fig.|5| we show the corresponding partial waves together with the mixing parameter ei for the best 
global fit. In some cases, the differences between NLO and NNLO are modest, in ^Pi and ^Pi NLO 
is even somewhat better. That means that the chiral TPEP is too strong in these phases. Note also 
that in the ^Pi phase OPEP is dominant. Thus, the inclusion of the contact interaction does not lead 
to a visible change. In ^P2, NNLO is still too strong but the prediction is considerably better than the 
NLO one. The energy dependence of ei is fairly precisely described at NLO and NNLO. These results 
are visibly better than the ones obtained in ref.[P] or in ref.|jl^, the latter being a NNLO calculation 
in the KSW scheme. This is shown in detail in fig.^, where ei is plotted versus the cms nucleon 
momentum {p < 350 MeV) in comparison to the Nimegen PSA and the results from refs.|12].^^ 



In the KSW approach, the pions are treated perturbatively. From these results we conclude that in 
the two-nucleon system, pions have to be treated non~perturbatively (if one intends to describe data 
above p ~ 150 MeV). 



*^^Note, however, that in order to fix the LECs in the ^5*1 — ^ Di channel we have used the leading coefficient gi = 1.66 
in the momentum expansion of the ei. This value agrees with the one given in ref. |28[ . In principle, we could take apart 
from the ^Si scattering length and effective range the deuteron binding energy as the third quantity to fix the three free 
parameters in the potential. However we refrain from doing that because of a strong correlation between the deuteron 
binding energy and the 5*1 effective range parameters. 

We are grateful to Iain Stewart for supplying us with the KSW NLO and NNLO results for larger momenta than 
given in refs.|ll|, |l2|. 
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5.1.3 D and F waves 

These are the most problematic phases in our approach. It is known that TPEP alone is too strong 
in some of these partial waves. To the order we are working, there are no contact interactions so that 
apart from varying the cut-off we have no freedom here. Still, for our best global fit, the D-waves 
are quite reasonably described, see fig-Q in particular ^Di and ^D^. The result for the latter phase 



is rather astonishing, since e.g. in the Bonn potential |30| the correlated two-pion exchange is very 
important in the description of ^D^. Such vrvr correlations are only implicitely present to the order we 
are working. They are hidden in the strengths of the LECs ci and C3 as explained in ref.Q. The very 
accurate (and parameter-free) description of the ^Di wave is of course important for the deuteron 
channel to be discussed in section [5^ . Even in ^D2, we get a fair description at NNLO. This partial 
wave is generally believed to be very sensitive to the A-resonance. We come back to this point when 
we discuss the NNLO-A approach. Note that the D~waves are rather sensitive to the choice of the 
regulator cut-off, see the upper panel in fig.^. At N^LO we have besides new pion-exchange terms 
the first contact interactions that feed into the D-waves, 

r^nlKont'^^r^n)-/^''. (5.4) 

More precisely, we have one independent parameter in each D-wave. Thus the cut-off independence 
will be restored by the running of the corresponding LECs. To illustrate this, we show in the lower 
panel of fig]^ the partial N'^LO results for the ^D2 channel. The corresponding potential consists of 
the NNLO terms plus one N^LO contact interaction. As expected, the cut-off dependence of the 
phase shift is very much reduced compared to the NNLO result. Of course, this illustrative example 
can not substitute for a complete N'^LO calculation, but one should expect very similar results. Note 
that according to our findings, the NNLO potential in the D-wave channels is not weak enough to be 
treated perturbatively, as it has been done in ref.[^. The potential has to be iterated to all orders in 
the LS equation. Only then one obtains a reasonable description of the phase shifts in these partial 
waves. Concerning the F-waves, which are shown in fig.y, ^F-^ and €3 are well described, whereas 
the NNLO TPEP is visibly too strong in ^^2, ^F^, and ^^4. This can be cured at higher orders by 
contact interactions. More precisely, a N'^LO calculation should be sufficient. 

5.1.4 Peripheral waves 

In figs.[To|, \n\ and ^ we show the G-, H- and I-waves together with the mixing parameters £4,5,6- 
These partial waves were first discussed in detail by the Munich group Q. Their calculation was 
perturbative and based on dimensional regularization of the TPE graphs. However, for these partial 
waves the iteration becomes unimportant and our findings confirm their results. The description of 
^Ga, ^Gs, 3^5, ^H^, ^He, 3/5 and ^Ij is visibly improved by the NNLO TPEP. Only in ^la the NLO 
result is better than the NNLO one. Of course, for the peripheral partial waves OPE does already 
a fairly good job, but the improvement in some of these phases due to the NNLO TPEP clearly 
underlines the importance of chiral symmetry in a precise description of low-energy nuclear physics. 

5.2 Deuteron properties 

We now turn to the bound state properties. At NNLO (NLO), we consider an exponential regulator 
with A = 1.05 (0.60) GeV, which reproduces the deuteron binding energy within an accuracy of about 
one third of a permille (2.5 percent). We make no attempt to reproduce this number with better 
precision. ^^ The results for the phase shifts, which correspond to these values of the exponential 



*^^Note that the deuteron binding energy is not used to fit the free parameters in the potential. 
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regulator, are very similar to those obtained with the sharp cutoff A = 
pleteness, we list in table |5| the values of the coupling constants in the '^S'l 
to the exponential regulator. 



0.875 (0.50) GeV. For com- 
-^Di channel corresponding 





C-iSi 


C^5i 


C-AD^-'J'Si 


NLO 

NNLO 

NNLO-A 


-0.0363 

-14.497 
-8.637 


0.186 
15.588 
7.264 


-0.190 
-4.358 

-0.447 



Table 5: The values of the LECs as determined from the ^Si —^ Di channel. We use an exponential 
cut-off with A = 0.6 GeV and 1.05 GeV at NLO and NNLO (NNLO-A), respectively. The LEG 
Cs^j is in 10^ GeV~^, while the others are in 10^ GeV~^. The parameters of the NLO, NNLO and 
NNLO-A potentials are obtained from fitting to the Nijmegen PSA. 



In table |6| we collect the deuteron properties in comparison to the data and two realistic potential 
model predictions (the pertinent formulae are given in app. |C|). We give the results for NLO and 
NNLO. We note that deviation of our prediction for the quadrupole moment compared to the empirical 
value slightly larger than for the realistic potentials. The asymptotic D/ S ratio, called r/, and the 
strength of the asymptotic wave function, As, are well described. The D-state probability, which is 
not an observable, is most sensitive to small variations in the cut-off. At NLO, it is comparable and at 
NNLO somewhat larger than obtained in the GD-Bonn or the Nijmegen-93 potential. This increased 
value of Pd is related to the strong NNLO TPEP. At N^LO, we expect this to be compensated 
by dimension four counterterms. Altogether, we find a much improved description of the deuteron 
as compared to ref.Q. Our results are almost as precise as the ones obtained in the much more 
complicated and less systematic meson-exchange models. 





NLO 


NNLO 


NNLO-A 


Nijm93 


GD-Bonn 


Exp. 


Ed [MeV] 


-2.1650 


-2.2238 


-2.1849 


-2.224575 


-2.224575 


-2.224575(9) 


Qd [fm2] 


0.266 


0.262 


0.268 


0.271 


0.270 


0.2859(3) 


T] 


0.0248 


0.0245 


0.0247 


0.0252 


0.0255 


0.0256(4) 


Td [fm] 


1.975 


1.967 


1.970 


1.968 


1.966 


1.9671(6) 


As [fm-i/2] 


0.866 


0.884 


0.873 


0.8845 


0.8845 


0.8846(16) 


Pd[%] 


3.62 


6.11 


5.00 


5.76 


4.83 


- 



Table 6: Deuteron properties derived from our chiral potential compared to two "realistic" potentials 
(Nijmegen~93 and GD-Bonn) and the data. Here, r^ is the root-mean-square matter radius. An 
exponential regulator with A = 600 MeV and A = 1.05 GeV at NLO and NNLO (NNLO-A), in order, 
is used. 



The coordinate space S- and D-state wave functions obtained in our approach are shown in fig.p!3|. 
At NLO they look quite similar to the ones obtained from various potential models. At NNLO one 
obtains a lot of structure in the wave functions below 2 fm. This is because two additional spurious 
(unphysical) very deeply bound states appear in the ^Si —^ Di channel. The binding energy of these 
states varies strongly by changing the cut-off. For the exponential regulator with A = 1.05 GeV we 
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get binding energies of Ei = 47.1 GeV and E2 = 2.5 GeV, respectively. This values correspond to 
center-of-mass momenta of about a few GeV which is clearly out of the applicability range of the 
low-momentum effective theory. Further, because of such huge values of the binding energy these 
unphysical states obviously do not influence physics in the energy region below 350 MeV that we are 
interested in and can, in principle, be integrated out.^^ A very similar situation with the spurious 



states happens at NNLO also in other S-, P- and D-waves. Note that these unphysical states are 



purely short range effects: as one can see from the fig. 14 the wave-functions corresponding to such 
states become negligeble for distances above 2 fm. The corresponding root-mean-square matter radii 
are (r^)^'^ = 0.27 fm and (r^)^'^ = 0.40 fm. For separations above 2 fm the NNLO deuteron wave- 



function is very close to the one obtained with the CD-Bonn potential. ^_ This is shown in figJlq. For 



a discussion on such deeply bound states in effective field theories, see ref.|23|. We end this paragraph 
with the following remark: According to Levinson's theorem, the difference between the phase shift 
at the origin and at infinity is given by nvr, with n the number of bound states. Thus, the phase 
shifts in the S-, P- and D-waves should become unphysical at large energies. This is, however, of 
no relevance for the EFT since we do not attempt to correctly reproduce (or predict) the phase shift 
behaviour for all energies (from threshold to infinity). 

5.3 Results for the NNLO-A approach 

As already stated before, our NNLO-A approach is not complete since we are omitting the dimension 
two pion-nucleon interactions not generated by the isobar. While the LECs C3 and C4 are dominated 
by the A, see ref. [Q], most of the correlated two-pion exchange is parametrized in ci. However, we are 
mostly interested in investigating the role of the A in all partial waves and therefore keep the cut-off 
A fixed at 875 MeV but refit the LECs Cj (see table |l] for their values). More precisely, a best global 
fit leads to a very similar cut-off value as in NNLO (this is expected due to the cut-off sensitivity of 
the D-waves) but for better comparison we discuss here the results obtained with exactly the same 
value for A. Note also that the precision of the fits is better than NLO in the theory without A 
but somewhat worse than the corresponding NNLO fits. This is due to the absence of contributions 
from higher mass states encoded in the LECs 01^3^4 not present in the NNLO-A approach discussed 
here. It goes beyond the scope of this paper to systematically include the effects of this important 



resonance in the framework of the EFT expansion as detailed in ref.|31|. In fact, a study of pion- 
nucleon scattering in that framework is not yet available and thus the corresponding LECs are not 
determined. Formally, we follow the Munich group [^2[ (for details, see app.y). Again, it is important 
to stress that we iterate our potential to all orders. 

Let us now discuss the results of the NNLO-A approach. We refrain from showing all partial waves 
but rather discuss some particular examples, collected in fig.|l^. The two S-waves shown in that figure 
are not very different from the NNLO result, although the description of "^Si is slighly worse at higher 
energies. All P-waves are very similar in NNLO and NNLO-A, the most visible difference appears in 
ei, as can be seen in the figure. The most dramatic effects appear in the D-waves. This is expected 
since these are parameter-free predictions and we had already pointed out the cut-off sensitivity in 
sec. ^.1.3 . Interestingly, the description of ^D2 is almost identical in the two approaches, consequently 



any important isobar effect in this partial wave can be well represented by contact interactions with 



*^^As soon as one is dealing with only the two-nucleon system there is no need to integrate out such unphysical states, 
since this does not modify the low-energy observables. We refrain here from the discussion of the complications which 
may arise in three- and more-body calculations due to such spurious states. Note, however, that according to our power 
counting one has an additional contact three-body force at NNLO which possibly can compensate the effects of the 
spurious states. 

*^ We would like to thank Hiroyuki Kamada for supplying us with the deuteron wave-function calculated with the 
CD-Bonn potential. 
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their strength given by the coupUng of the A to the irN system. In "^D^ (also sliown in fig. 16) 
the absence of the scalar-isoscalar two-pion correlations is clearly visible. Our result thus confirms 
a finding made already in the Bonn potential, namely that this particular partial wave is essential 
dominated by correlated TPE. We also note that the description of ^Di and ^D2 is worse in NNLO-A 
than in NNLO. For the F-waves, the differences between the two approaches are very small, with the 
exception of ^F^, which is improved in the presence of the isobar. The most significant effect in the 
higher partial waves shows up in ^G^ as shown in fig.|lB. Clearly, two-pion correlations not related 
to the A play an important role to bring the prediction for this partial wave in agreement with the 
data. The deuteron properties for an exponential regulator with A = 1.05 GeV are listed in table p. 
Most observables come out in between the NLO and NNLO results for the deltaless approach. The 
sole exception is the quadrupole moment, which is improved. We remark again that only after a fully 
systematic inclusion of the A, one should draw quantitative conclusions. It appears, however, that 
simply using resonance saturation to encode the phyiscs of the isobar in the dimension two pion- 
nulceon LECs is a sufficient procedure in the two-nucleon sector. A truely quantitative inclusion 
of the delta should be done in the framework of the "small scale expansion" [^l|. In that case, the 
leading delta effects would already appear at NLO since the nucleon-delta mass splitting is counted 
as a small parameter. Such an investigation is, however, beyond the scope of this article. We also 
note that explicit isobar degrees of freedom were considered in ref. [H] . A direct comparison with their 
work is, however, not possible since no separate result for fits with and without delta are given. 

6 Summary 

In this paper, we have calculated properties of the two-nucleon system based on a chiral effective 
field theory. The underlying formalism was already presented in I. The results of this investigation 
can be summarized as follows: 

1) Based on a modified Weinberg power counting (as explained in I), we have constructed a chiral 
two-nucleon potential at NNLO. It consists of one- and two-pion exchange diagrams, including 
dimension two insertions from the effective pion-nucleon Lagrangian. The corresponding LECs 
have been taken from an investigation of ttN scattering Q- In addition, there are two and seven 
four-nucleon contact interactions at LO and NLO, respectively. The coupling constants of these 
terms must be fixed by a fit to data. 

2) For large momenta, the potential becomes unphysical and has to be regularized. We perform 
this regularization on the level of the Lippmann-Schwinger equation, as explained in sectj2| 
using either a sharp or an exponential regulator function. At NLO, physics does not depend on 
the cut-off in the range between 400 and 650 MeV. At NNLO, this range is larger and extends 
from 650 to 1000 MeV. This can be understood from the chiral TPEP, which at NNLO includes 
vrvr correlations. These introduce a new mass scale well above twice the pion mass. 

3) We have shown that the contact interactions can be combined in such a way that each com- 
bination feeds into one partial wave. More precisely, the nine four-nucleon couplings can be 
determined uniquely by a fit to the two S-waves, four P-waves and the mixing parameter ei for 
nucleon laboratory energies below 100 MeV. As expected from the power counting underlying 
the EFT, the fits improve when going from LO to NLO to NNLO, compare fig.|l[ 

4) At NNLO, the resulting S-waves are of very high precision (for nucleon laboratory energies below 
300 MeV), see e.g. tables |2|,^ and fig.§ The so-called range parameters collected in table |^ 
agree with what is found in the phase shift analysis. The P-waves are mostly well described, in 
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particular the mixing parameter ei is in good agreement with the phase shift analysis. We also 
note that above nucleon cms momenta of about 150 MeV, our NLO and NNLO results are far 
better than the one obtained in the KSW scheme at NLO and NNLO. 

5) All other partial waves are free of parameters. The D-waves, in particular ^Di and ^D^ are 
very well described. We have also discussed the cut-off sensitivity of these results. The NNLO 
TPEP is too strong in the triplet F-waves. This is expected to be cured at N^LO due to the 
appearance of dimension four contact terms. For the peripheral waves, we recover the results 
of the Munich group P], namely that in most cases OPE works well but chiral NNLO TPEP 

h. 



clearly improves the description of some partial waves like e.g. ^Gs, ^H^ or ^/y 



6) The deuteron properties are mostly well described, at NLO and NNLO, compare table ^. At 
NNLO, the deuteron wave functions shows some interesting structure due to the appearance of 
two very deeply bound states. These are an artefact of the NNLO approximation. They have 
no influence on low energy properties and can be completely projected out from the theory. Our 
precise deuteron wavefunctions can be used for pion photoproduction, pion-deuteron scattering 



or Compton scattering off deuterium (still, the hybrid approach proposed by Weinberg [33| 
remains a useful tool). 

7) We have also considered an approach with explicit A degrees of freedom in the TPEP. This 
NNLO-A approach is very similar to the NNLO results in the theory without isobars, with the 
exception of the partial waves that are sensitive to pionic scalar-isoscalar correlations like e.g. 
^D'^. We conclude that the inclusion of the A via resonance saturation of the dimension two ttN 
LECs captures the essential physics of the isobar in the two-nucleon system. We note, however, 
that a more systematic study of pion-nucleon scattering in an EFT including the A is needed 
to further quantify these statements. 

Our findings do not only show that the scheme originally proposed by Weinberg works quantitatively, 
it even works much better than it was expected. It extends the succesfull applications of effective 
field theory (chiral perturbation theory) in the pion and pion-nucleon sectors to systems with more 
than one nucleon. Clearly, one should now reconsider processes, which have been evaluated using 



Weinberg's hybrid approach |3^] (vr — d scattering [33, ^], 7^ — > t: d [35|, 7^ — > 7^ ^^) and extend 
these considerations to systems with more than two nucleons. First steps in that direction within the 
potential approach can be found in recent works [R^, p^]. In addition, a fresh look at charge symmetry 



and charge independence breaking is called for (for earlier studies, see e.g. refs.|39, 4C, |41[|). Work 
along these lines is underway. Furthermore, the precise relation to the KSW scheme, which has been 
shown to be successfull at low energies, has to be worked out. 
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A The renormalization procedure 

In this appendix we spell out some details about the renormalization procedure mentioned briefly in 
sect.^. We will use the divergent loop integrals Jo,2,4 defined in eq.( |2.lD to express all divergences and 
remove these by proper subtractions and redefinitions of physical quantities. 
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We consider first the one-loop corrections to the OPEP. These were given in eq.(4.33) of I. We rewrite 
this expression here using different conventions for the isospin operators and the pion decay constant 

+2f {5i ■ q) {a2 ■ q)\ , 



(2/^)4 3 UJ^q J (27r)'^ Lof 

4^^ (fi • f2) (ai • q) {&2 ■ q) ^ JSM^ + 3M2 In ^ + 4J2 - QMIJ^ 



9a 



in terms of the two divergent loop functions Jo,2 defined in eq. (|2.1| ). Therefore, this contribution has 
exactly the same form as the OPEP (renormalized OPE), 

VS^ + VS> = Vl^ , (A.2) 

provided we redefine the coupling constant g\ (the superscript "0" denotes the leading term in the 
chiral expansion) in the following way: 



{9\? = {9\? ^'° '' 



^2 I 5M2 + 3M2 In ^ + 4 J2 - 6Af 2 Jo | . (A.3) 



12vr2/^ 

Clearly, g'j^ and 5^ differ by terms of second order in the chiral dimension. Consequently, all NLO 
one-loop corrections to OPEP can be taken care off by renormalization of g\. 

In addition, there are the one-loop corrections to the lowest order four-fermion interactions, cf. 
eq.(4.34) of I. Performing spin averaging, the corresponding contribution can be expressed in terms 
of the divergent loop functions Jo,2 as 

.,2 r 

^l^'ll'i:, ,-3 



(P ( JVf2 

= -^^^ C^T (3 - n • T2) (ai • as) JSM^ + ZMI In ^ + 4J2 - QM^h 



(A.4) 



or in a more compact notation 

'NN,l-loop 



^nS.i-Ioop = -35 ('^i ■^2) + S {a, ■ a2)(Ti • r^) , (A.5) 



with 

,2 



S = ^^ C^T {5M2 + 3M2 In ^ + 4J2 - 6M2 Joj . (A.6) 



Antisymmetrization allows to map the two spin-isospin operators appearing in eq.( A.5 ) onto the two 



non-derivative operators used in Vf^j^, cf. eq.(^^). Therefore, the effect of the one-loop corrections 



to the lowest order four-nucleon contact interactions can be completely absorbed by renormalizing 
the constants C^ and Cj<, 

CJ = C70 - 35 , C^ = C70 - 35 . (A.7) 
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We note that further renormalization of these couplings is due to the two-pion exchange, as dicussed 
below. 

A similar procedure also works for the one-loop TPE graphs, cf. eq.(4.30) of I. In fact, the NLO 
TPEP renormalizes the various dimension zero and two coupling constants Cj. So as not to repeat 
the argument, we simply give the relevant momentum space integrals in terms of J0.2 and J4, 



1 



dH 



(27r)^ a;+a;_ (a;_(_ + w_ ) 

dH P 

(27r)'^ u;_l_a;_ (a;_(_ + w_ ) 



dH l^ + q^ 



dH {q • If 

(27r)'^ cj-l-tj- (a;_|- + tj_ ) 

dH {UJ+ - UJ^f 

(27r)'^ u;_l_a;_ (a;_|_ + w_ ) 



87r2 \ 9 s- 



M2 

In ^ - 2 Jo 



(A. 



1 



247r2 



2s3 

m 

q s-q 



s + Q . ...,.2 „2^ , o^n«^2 , „2^^ M^ 



+ (14M^ - q') + 2(9M^ + q') In ■ 



+6J2-4(9M2+g2)jA 



(A. 



^ / — (512M6 + 384M^g2 ^ isgM^/ + 235^) In ^^ 



247r2 [5gs 
+ — (-898M^ + 192M2g2 + ^4^ _ -(450M^ + SOM^g^ + 23g^) In 



Mi 



+ 6J4 - 4(9M2 + g2)j2 + ^(450M^ + 50M2g2 + 23g^) Joj , 



-^ 2My In ^ - ^ + -g^M^ + -g2 J2 - 4M2g2 Jo 



87r2 

1 
6^ 



— In 



s + 



5 3 

n M^ 9 9 

g2ln^ + 8M2 + 2g2jo 



(A.IO) 

(A.ll) 
(A.12) 



with s = y^AM'^ + g2 and 1^-1- given in eq.(4.24) of I. All other integrals appearing in the 1-loop NLO 
TPEP can be deduced from these expressions by taking proper linear combinations or differentiation 
with respect to M^. Putting pieces together, the expression for ^"2^ ^_ioop takes the form 



V. 



(2) 



27r,l— loop 



taTPEP 
•'NLO 



+ {Si + S2 q^) (ri • fa) + ^3 



((Ti • q) ((T2 • q) - (ai • (T2) q' 



(A.13) 



with 



^2 
^3 



(A.14) 



^^l^^^-lSMl{bg\ - 2g\) In ^ - Ml{%lg\ - Ug\ + 4) 

+iml{bg\ - 2g\)Jo - 3{3g\ - 2g\)J2 
384^2/4 {(-23gA + 105a + 1) In ^ - {I2,g\ + 2g\) + {22,g\ - lQg\ - 1) Joj, (A.15) 



and the non-polynomial part I^^lo"^ is given in eq.(2.12). The polynomial terms clearly renormalize 
the coupling constants of the dimension zero and two four-nucleon contact terms. Again, we perform 
antisymmetrization to map the terms appearing in eq.( A.13| ) onto the basis used in V^;^ , cf. eq. (^^) . 
Including the contribution from eq.(A.7), the complete renormalization of the four-nucleon couplings 
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takes the form 

/^^r /^^O c /^^^ /^^ A Q /^^ /^O C 

Oi — Oi — O2 , U2 — '-'2 ~ ^'-'2 1 '-'3 — ^3 ~ '-'3 1 

CI = cl - 452 , CI = c'i , Q" = cl + S3 , a; = c^j . (a.i?) 

Note that C5 and C7 do not get renormahzed to this order. All coupling constants appearing in the 
main text are understood as the renormahzed quantities discussed here. 



B Inclusion of the A(1232) 



The inclusion of the A(1232) in the TPEP has been worked out in ref.[32|. For completeness, we 
collect here the pertinent formulae. There are three distint contributions. 



A-excitation in the triangle graphs: 



with 



aTPEP _ 

A, triangle 

L{q) 


-I^(n-2){(6i?-^V(?) + 12A^ 
s s + q 


D{q) 

A 

E 


A J2M^ IJ? + q^ 2A 
= mA - m = 293 MeV , 

= 2Ml + q^ - 2A2 . 



Single A-excitation in the box graphs: 



^^ (n • T2) ((12A2 - 2QMl - llq^)L{q) + 6E^D{q) 



I92n^n 



with A{q) given in eq.(2.13|). 



Double A-excitation in the box graphs: 

^4 



^Z^bi^-2 = --^^[-4A'L{q) + E[H{q)+{E + 8A')Diq)]^ 



^^ (n • T2) 1{12E - s^)L{q) + 3E[H{q) + (SA^ - E)D{q)] 



384^2;^ 



(B.l) 



(B.2) 

(B.3) 

(B.4) 
(B.5) 



128^2/4 ((^1 • 9"')(^2 • q) - q\o, ■ 0^2)) {-2L(g) + {s^ - AA^)D{q)^ 
Y^^I^ (n • r2) ((ai • q){a2 • q) - q\a, • ^2)) s' A{q) , (B.6) 



J^^TJI (^1 • ^2) (('?! • q){ff2 ■ q) - q\S, ■ a^)) [2L{q) + {AA? + s^)D{q)^ , 
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with 



H{q) 



2E 



4A^ 



L(g)-L(2v/A2-M2) 



(B. 



Finally, let us show how to arrive at the coupling constant relation given in eq.( |2.16| ). In spin-isopsin 
SU(4) (or more generally, SU(6)), one obtains 



SAAfTT = -j= g-KN 



Using now the Goldberger-Treiman relation 



m 
OttN = gA^r 



(B.9) 



(B.IO) 



one arrives at eq.(2.16). The numerical value of ^aa^tt follows from using g^ = 1.26. Note that this 
value for the coupling constant leads to a A-width of about 117 MeV. 



C Formulae for the deuteron properties 

Here, we collect the formulae needed to calculate the deuteron properties. Denote by u{r) and w{r) 
the S- and D-wave coordinate space wave functions. We denote the momentum space representations 
of u{r)/r and w{r)/r by u{p) and w)(p), respectively. We have 



Normalization 
D — state probability 
Quadrupole moment 



Asymptotic S — state 
Asymptotic D/S — ratio tj 







dpp ['u{p) + 'w{p) ] = dr [u{r) + w{r) ] = 1 

Jo 



dpp w{p) = / drw{r) 
Jo 

1 f°° 
Qd=7^ drr'^w{r)[V8u{r) 

^U Jo 



w(r)] 



YoL 'n^p dp dp +'^"(^^ 

„2 



du{p) 
dp 



(^^\ + euW 



^'\:^) 



u{r) 



Ase ^"^ for r 



w[r) 



RMS (matter) radius : r^ 



oo , 

7?A5(1 + — + — -2 

7r (7rj^ 



drr [u{r) + w{ry 



-jr 



1/2 



for r 



oo 



(C.l) 
(C.2) 



(C.3) 
(C.4) 
(C.5) 

(C.6) 



with 7 = y/m \Ed\ = 45.7 MeV (using m = (mp+m„)/2). Note that the momentum-space representa- 
tion of Qd given in eq.(C.3) shows why one cannot use a sharp momentum-space regulator to calculate 
this quantity. We also remark that the D-state probability is not an observable. Meson-exchange 
current corrections to Qd are not given. 
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Figures 
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Figure 1: Quadratic deviations in the phase shifts, {6 — 5 )^, versus iiJiab in GeV. Data from the 
Nijmegen PSA have been fitted for ii^iab < 0.1 GeV. The purple dotted, green dashed and red sohd 
curves represent LO, NLO and NNLO results, respectively. 
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Figure 2: Running of the four-nucleon coupling constants in the ^Sq (upper 
panel) and ^Si —^Di (lower panel) partial waves. In the upper panel, the solid 
red (dashed blue) line refers to Cig^^ (Cis^). In the lower panel, the solid red 
(dashed blue) line refers to Cag^ (Cas^)- The green (dashed-dotted) line refers 
to the constant in the coupled ^Si — ^ Di system. The units are the same as 
in the table ||. 
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Figure 3: Phase shifts for the two solutions in the ^5o-wave as dicussed in 
the text. At NLO (upper panel), these are indistinguishable. At NNLO (lower 
panel), one of the solutions (black line) shows an unacceptable behaviour at 
higher momenta and is discarded. 
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Figure 4: Predictions for the S-waves (in degrees) for nucleon laboratory 
energies E'lab befow 300 MeV (0.3 GeV). The purple dotted, green dashed and 
red solid curves represent LO, NLO and NNLO results, in order. The blue 
squares depict the Nijmegen PSA results. In the upper and the lower panel, 
^Sq and ^Si, respectively, are shown. 
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Figure 5: Predictions for the P-waves and the mixing parameter ei (in degrees) for nucleon laboratory 
energies E'lab below 300 MeV. For notations, see fig.0. 
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Figure 6: Predictions for the mixing parameter ei for nucleon cms momenta p befow 
350 MeV. The green, blue and red sohd curves represent our LO, NLO and NNLO results, 
in order. For comparison, the NLO O] and NNLO |12] results in the KSW scheme are 
also shown. The black squares depict the Nijmegen PSA results. 
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Figure 7: Predictions for the D-waves and the mixing parameter e2 (in degrees) for nucleon laboratory 
energies E'lab below 300 MeV. For notations, see fig.0. 
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Figure 8: The ^D2 partial wave at NNLO for three different values of the 
cut-off A (upper panel) and at N^LO (lower panel). The blue diamonds are 
the result of the Nijmegen partial wave analysis. 
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Figure 9: Predictions for the F~waves and the mixing parameter £3 (in degrees) for nucleon laboratory 
energies E'lab below 300 MeV. For notations, see fig.0. 
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Figure 10: Predictions for the G-waves and the mixing parameter £4 (in degrees) for nucleon 
laboratory energies -Eiab below 300 MeV. For notations, see figj^. 
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Figure 11: Predictions for the H-waves and the mixing parameter es (in degrees) for nucleon 
laboratory energies £^iab below 300 MeV. For notations, see figj^. 
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Figure 12: Predictions for the I-waves and the mixing parameter ee (in degrees) for nucleon labora- 
tory energies -Eiab below 300 MeV. For notations, see figJ3. 
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Figure 13: Coordinate space representations of the S- (red solid line) and 
D-wave (blue dashed hne) deuteron wave functions at NLO (upper panel) and 
NNLO (lower panel). 
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Figure 14: Coordinate space representations of the S~ (red solid line) and 
D-wave (blue dashed hne) for the two unphysical boundstates. 
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Figure 15: Coordinate space representations of the S- (upper panel) and 
D-wave (lower panel) deuteron wave functions at NNLO compared to the one 
from the CD-Bonn potential. 
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Figure 16: Phase shifts (in degrees) for the NNLO-A theory (red sohd hnes) in comparison to the 
NNLO (green dashed hnes) and NLO (purple dotted hnes) resuhs as a function of the lab energy in 
GeV. For other notations, see fig]^. 
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